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Abstract 

We study the substrate inhomogeneity effect on the wetting of a spin- 1/2 Ising 
ferromagnetic film in an external magnetic field H, using Monte Carlo simulations. It is 
found that the inhomogeneity leads to the formation of clusters of positive and negative 
spins in each layer. However, depending on the values of H, for a fixed surface magnetic 
field H s , each layer exhibits three different phases; totally wet, no wet and partially 
wet. In the latter case, we show the existence of three distinct configurations namely: A 
configuration in which the layer is partially wet with a total disconnection between clusters 
(PWTD); a configuration in which the layer is partially wet with a partial disconnection 
between clusters (PWPD); the configuration in which the layer is partially wet and the 
clusters are totally connected (PWTC). Furthermore, we show that an increase of H s 
values breaks the bonds connecting some clusters of the phase (PWTC) which leads to 
an expansion of the (PWPD) region. On the other hand, the frequency distribution and 
the average size of clusters are investigated in the (PWTD) region for fixed values of 
temperature T, H s and H. Moreover, we show the existence of 3D-islands resulting from 
the formation of islands in each layer. 
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1 Introduction 



Multilayer films adsorbed on attractive substrates may exhibit a variety of possible phase 
transitions, as has been reviewed by Pandit et al. [1], Pandit and Wortis [2], Ebner et al. [3] 
and Patrykiejew et al. [4]. The existence of phase diagrams depending on the strength of the 
substrate potential in a lattice gas model for multilayer adsorption, have been shown by Beny- 
oussef and Ez-Zahraouy [5]. The wetting and layering transitions are caused by the competition 
between short-ranged attractions and long-ranged repulsions in systems like langmuir mono- 
layers and adsorbed monolayers [6] . These transitions have also been found, under the effect of 
the geometry of the system, in the case of the double wedge wetting [7] in magnetic films, and 
the wetting behavior of associating binary mixtures at attractive walls [8] . We showed recently, 
in Ref. [9], the existence of these transitions under the effect of the surface coupling, when 
using the mean field theory. 

On the other hand, many recent studies showed that indirect interactions can lead to the for- 
mation of nano-structures at surfaces [10] mediated by the structure and the geometry of the 
film. Monte Carlo Simulations [11] and subsequently derived mean field theories [12], showed 
the existence of higher island densities than those expected by standard nucleation theory. It 
is shown that several experimental works, e.g. [13], supports these predictions, such as the 
deposition of Au on mica substrates at a high temperature [14]. 

A small variation of the preparation conditions may considerably change the obtained magnetic 
structure and allows the investigation of atomic structure as well as the magnetic domain struc- 
ture of several nanostructured systems [15,16]. However, Lenz and Lipowski [17] showed that 
the morphology of wetting layers on structured surfaces is determined by the geometry of the 
underlying surface domains, and that the wetting layer exhibits several distinct morphologies. 
On the other hand, it has been shown [18] in the case of a two-dimensional lattice-gas model, 
that some basic regimes of drying mediated nanoparticle assembly are distinguished by the 
spatial uniformity of solvent dynamics, and by the fluctuations of nanoparticle domain 



2 



boundaries following evaporation. It has also shown that when solvent disappears 
homogeneously from the surface, disk-like or ribbon- like domains reminiscent of spinoidal 
decomposition form at early time. If instead domain boundaries are frozen following 
evaporation, dynamical constraints arrest this growth at an early stage. Network structures 
are formed as vapour nuclei meet [18], otherwise networks fragment to form distinct 
domain that asymmetrically evolve as in homogeneous coarsening. 

Our aim in this work is to study the effect of the inhomogeneity of substrate on the wetting 
transitions of a spin- 1/2 Ising film when a surface magnetic field is applied on alternate clusters 
of the surface. The paper is organized as follows. In section 2, we describe the model and the 
method used: Monte Carlo (MC) simulations. In section 3 we present results and discussions. 

2 Model and Monte Carlo simulations 

Our model can describe either a magnetic system in which the Ising spin variables are 
Si = +1,-1, or a physical state of the liquid-gas model where the occupation variables 

rii = +1, (+1 when the site is occupied and if the site is not occupied) can replace 
the spin variables when rii = (1 + Si)/2. The external magnetic field H and the surface 
magnetic field H s i can then play the role of the chemical potential and the substrate potential, 

respectively. These analogies concerning the lattice gaz model and a magnetic system are more 
explained in Refs. [1-4]. The case a two-dimensional lattice gas model is studied in details in [18]. 
The system we are studying, illustrated by Fig. 1, is a magnetic thin film formed with N = 4 
layers coupled ferromagnetically. Each layer is a square of dimension N x x N y = 64 x 64 spins. 
N x and N y stand for the number of spins in the x and y directions, respectively. A surface 
magnetic field H s is acting only on alternate clusters of spins of the surface k — 1. This is 
represented, in Fig. 1, by symbols (+) for alternate clusters of dimension l p x l p = 8 x 8. 
Whereas H s is absent for the remaining alternate clusters of dimension l n x l n = 8 x 8. These 
clusters are represented by symbols (o). 
The Hamiltonian governing this system is given by 

w = - E JijSiSj - E(# + H ai )Si (i) 

<i,j> i 
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where, Si(l = = — 1, +1 are the random spin variables and 
the interaction between different spins is assumed to be constant Jij = J. 
The first summation in Eq. (1) includes only the nearest neighbor spins. The second 
summation runs over all sites i of the system. H is an external magnetic field applied 
on each site of the film. 

The surface magnetic field Hsi, appearing in the Hamiltonian (1), is applied only on sites % of 
the surface k — 1, and distributed alternatively, according to: 



The inhomogeneity in the surface results from this distribution of the magnetic field Hsi. 
We use the Monte Carlo simulations under the Metropolis algorithm. 
An attempt to flip a randomly chosen spin in site '%' from up to down, or vice versa, 
is envisaged using the Metropolis probability: 



where, AE is the resulting change in energy after the flip of the spin V, ks is 
the Boltzmann constant and T the absolute temperature. The flip of the spin is accepted 
if p M is greater than a random number p r generated by the computer so that p r e[0, 1]. 
Otherwise the flip of the spin located at the site 'i' is ignored if the randomly generated 
number is so that p r > pu- One Monte Carlo step (MCS) is reached once the number 
of visited sites is close to the number of sites of the system. For a fixed temperature and/or 
surface magnetic field, runs of 5 x 10 6 Monte Carlo steps (MCS) are performed for 
different initial conditions. Mean values of computed quantities are estimated for different 
initial conditions, with the discard of the first 5 x 10 5 (MCS) for each run. 
A preliminary study showed that the topologies of the phase diagrams are similar (with dif- 
ferent numerical values) when the film thickness varies from N = 3 to N = 9 layers, when 
the number of spins of each layer varies from N x — N y — 16 to N x — N y — 64; where 
N x and N y are the number, of spins of each layer, in the x and y— directions, respectively. 
When increasing the film thickness, the topology of the established phase diagrams are 
similar for thin films (N < 10) layers. For simplicity and fast Monte Carlo simulations, and 
taking into account the above considerations we give numerical results in this paper for a film 




+H S for all sites i e clusters with symbols (+) 
for all sites i e clusters with symbols (o) . 



(2) 



p M = min[l,exp(—AE/k B T)], 



(3) 
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with N = A layers and N x — N y — 64 spins for each layer. 

We mean by an 'island', a set of negative spins surrounded by positive ones. In all the following, 
our interest will be concentrated on the mean size of islands (MSI) and the average number of 
islands (ANI). The different phases are determined by values of the average number of islands 
and their mean sizes as following: 

The PWTD phase is characterized by: ANI < 32 and MSI < (MSI < N x x N y )/2. 
The PWPD phase is present for: ANI < 32 and (N x x N y )/2 < MSI < N x x N y . 
The PWTC phase corresponds to: ANI > 1 and MSI > (N x x N y )/2. 

However, in the particular case of the non wet phase, we have ANI = 1 and MSI = Nx x Ny 
(which is the total number of spins of each layer). We note that in this particular case , although 
positive spins are completely absent in the layer, we assume that this situation corresponds to 
a unique island. 

3 Results and discussion 

A sketch of the geometry of the system we are studying, is presented in Fig. 1. In order to 
outline the behaviour of the mean size of islands (MSI) and the average number of islands (ANI) 
as a function of the external magnetic field, we plot in Figs. 2a and 2b, these parameters, for 
a surface magnetic field H s = 2.5 and a temperature T = 3.5. Indeed, Fig. 2a corresponds 
to the surface k — 1, whereas Fig. 2b is plotted for the last layer k = N = 4. It is found 
that, for the surface k — 1 (Fig. 2a) and the last layer k = N (Fig. 2b), when decreasing the 
external magnetic field from H = 0, the average number of islands (ANI) exhibits two islands 
for H ~ —0.1, with a rapid increasing of the number of negative spins in each layer, as it is 
illustrated by the mean size of islands (MSI) curve. For T = 3.5 and —0.16 < H < —0.07, 
the average number of islands (ANI) increases rapidly and the mean size of islands (MSI) un- 
dergoes a local maximum. When decreasing the magnetic field H so that H < —0.20 the 
average number of islands of the surface stabilizes at 32, whereas it drops suddenly to one 
cluster formed with negative spins, in the case of the last layer k = N. In the case of the 
surface, the mean size of islands (MSI) undergoes a large local plateau corresponding to the 
value 64, followed by the apparition of a unique island for H ks —0.80. This scenario is ab- 
sent in the case of the last layer k = N, Fig. 2b, where the islands are rapidly connected to 
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each other for H < —0.15. On the other hand, the increasing surface magnetic field effect, on 
island formation is illustrated by Figs. 3a and 3b for a fixed temperature T = 3.5 for several 
values of the surface magnetic field: H s = 1.5, H s = 2.0 and H s = 2.5. Indeed, for small 
values of the external magnetic field, —0.20 < H < 0, the increasing surface magnetic field 
H s is not felt at this "higher" temperature value 3.5. While, for H < —0.20, the increasing 
surface magnetic field amplitude is to produce a large step of mean size island value 64, see 
Fig. 3a. It is worth to note that this step disappears for very small surface magnetic field 
values. Concerning the average number islands (ANI), the same arguments still valid when 
increasing the surface magnetic field amplitude H S) Fig. 3b. The islands are rapidly connected 
to each other for very small values of H s . On the other hand, we show in Fig. 4 the existence 
of three regions of island configurations in the plane (T,H), for k = N and a fixed surface 
magnetic field H s = 2.0. These three different regions are: the totally wet (TW) region, the no 
wet region (NW) and the partially wet region. In the latter case, see in Fig. 4, the existence 
of three distinct island configurations, namely; the (PWTD) configurations where the layer is 
partially wet with a total disconnection of the islands; the (PWPD) configurations in which the 
layer is partially wet with a partial disconnection of the islands; The (PWTC) configurations 
corresponding to a partially wet layer with totally connected clusters. Indeed, at low temper- 
ature values, only the totally wet (TW) phase for H 0, and the non wet (NW) phase for 
H < —0.5 are present. At higher temperature values, provided that T is kept less than the 
critical value T c = 3.87, the partially wet phase occurs. Starting from the totally wet phase 
(TW), present for positive values of H, the decreasing of the external magnetic field (from to 
negative values), gives rise to the phase (PWTD) showing totally disconnected configuration 
of islands. This phase is followed by the (PWPD) phase in which the configuration of islands 
are partially disconnected when —0.16 < H < —0.07. For H < —0.20, the partially wet phase 
(PWTC) where the islands are totally connected to each other, is reached. Decreasing the 
magnetic field H more and more, the non wet phase (NW) is found and persists for H — > — oo. 
The topology of the established phase diagrams stand valid for the other layers k < N 
of the film. The layer k=4 is chosen because the studied phases in the partially wetting region: 
totally disconnected (PWTD), partially disconnected (PWPD), and totally connected (PWTC) 
are well outlined in this layer. Moreover this layer is chosen in order to outline the range effect 
of the substrate on deeper layers. 
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In order to more clarify these findings, we plot in Figs. 5a and 5b, different island configurations 
in different regions of the phase diagram for the different layers of the film. 
Indeed, Fig. 5a, shows that, for H s = 2.0, T = 3.5 and H = —0.13, the first layer k = 1 
shows a configuration corresponding to a partially wet layer with totally connected clusters 
(PWTC). The second layer k = 2 belongs to the partially disconnected configuration (PWPD), 
whereas the third layer k = 3 is in the totally disconnected phase (PWTD). On the other hand, 
when maintaining the parameters T and H s constant, and varying only the external magnetic 
field, a given layer undergoes different configurations belonging to different regions: totally 
connected (PWTC), partially disconnected (PWPD) and totally disconnected (PWTD). This 
situation is illustrated by Fig. 5b for the last layer k = 4. Indeed, for H = —0.13 this 
layer is in the totally connected phase (PWTC). For H = —0.15 the configuration of this 
layer partially disconnected (PWPD); and when H = —0.17 the islands of this layers are com- 
pletely disconnected from each other so that this configuration belongs to totally disconnected 
(PWTD) region. The numerical values given in this work correspond to N x — N y — 64. 
However, a preliminary study showed that when N x — N y — 32 and T = 3.5 the partially 
disconnected region (PWPD) is located in the interval —0.13 < H < —0.09. Hence, the finite 
size effect and thermodynamic limit is to enlarge partially wetting regions (PW). Moreover, 
the formation of islands in each layer leads to the formation of the three dimensional islands 
(3L>-islands). The superposition of these 2D-island of adjacent layers leads to the formation of 
such 3-D-islands in the film. For example, Fig.5a shows these 3-D-islands, when superposing the 
layers k — 1, 2, 3 over each other. This situation is more outlined, for example, in the rectangle 
of coordinates in the rectangle 16 < x < 24 and 32 < y < 40 for the layers k = 1,2,3 and 4, re- 
spectively. The superposition of these islands constitute a three dimensional island in the bulk 
of the film. However, the height of the formed islands depends on the values of temperature, 
surface and external magnetic fields. In the example of Fig. 5a, some island heights are close to 
4. A preliminary study showed that the island height decreases with decreasing values of H s , 
for fixed values of temperature T and external magnetic field H. To complete this study, the 
frequency distribution of different configurations, are plotted in Fig. 6. This figure shows the 
corresponding frequency island size distributions, for the totally disconnected region (PWTD), 
for the layers k = l, k = 2, k = 3 and k = 4 when the temperature is fixed at T = 3.5, 
the surface magnetic field is H s = 2.0 and H = —0.142. The different scenarios presented by 
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the system concerning the inhomogeneity substrate effect on the distribution of cluster sizes 
and number of existing clusters are presented. Indeed, the mean size cluster of negative spins 
undergoes a local maximum before exhibiting first and second steps. The plateau of the first 
step increases with increasing surface magnetic field values. 

4 Conclusion 

We have studied the effect of the inhomogeneity of a substrate on the wetting transitions of 
a spin-1/2 ferromagnetic Ising thin film under the effect of an alternate surface magnetic field 
H s acting on alternate clusters of the surface, using Monte Carlo simulations. In the par- 
tial wetting region each layer exhibits three different configurations namely: partial wetting 
with totally disconnected cluster (PWTD), partial wetting with partially disconnected clusters 
(PWPD) and partial wetting with totally connected clusters (PWTC). On the other hand, we 
found that increasing the surface magnetic field leads to the disconnection of some connected 
clusters in the (PWPD) region. Furthermore the 3-D-islands occur for strong surface magnetic 
field values. Moreover, the distributions of the cluster size, in each layer of the film, are also 
computed. 
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Figure Captions 

Figure 1.: 

A geometry sketch of the studied system formed with N 
is acting on alternate cluster spins of the surface k — 1. 
'+', and absent elsewhere: islands with symbols V. 

Figure 2: 

The dependency of the mean size of islands (MSI), and the average number of islands (ANI) of 
the surface k — 1 (a), and the last layer k = N, as a function of the external magnetic field H 
for a fixed temperature T = 3.5 and a surface magnetic field value H s = 2.5. 



= 4 layers. A surface magnetic field H s 
H s is present on islands with symbols: 
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Figure 3: 

(a) The mean size of islands (MSI), and (b) the average number of islands (ANI) of the surface 
k = 1 as a function of H for T = 3.5 and several values of H s : 1.5; 2.0 2.5. 

Figure 4: 

The (T, H) phase diagram for the layer k = 4, as a function of the external magnetic field H, for 
several values of the surface magnetic field H s : 1.5; 2.0 and 2.5. The regions corresponding to 
the different configurations: (TW) totally wet, (PWTC) partially wet with totally connected 
clusters, (PWPD) partially wet with partially disconnected clusters, (PWTD) partially wet 
with totally disconnected clusters, and (NW) non wet, are presented. 

Figure 5: 

island maps of positive spins for different layers of the film with N = A a) The first layer k = 1 
is in the (PWTC) configuration, the second layer k = 2 is in the configuration (PWPD), while 
the layer k = 3 belongs to the configuration (PWTD); for T = 3.5, H = -0.13 and H s = 2.0. 
b)For H = —0.13, H = —0.15 and H = —0.17, the last layer k = 4 exhibits three different 
configurations (PWTC), (PWTD) and (PWTD), respectively, for T = 3.7 and H s = 2.0. 

Figure 6: 

Frequency island size distributions for the layers k — 1, k — 2, k — 3 and k = 4 for the 
configuration (PWTD) corresponding to T = 3.5, H s = 2.0 and H = —0.142. 
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Fig. 1 




Hs=+H s ( symbols : + ) 
Hs= ( symbols : o ) 




H 



Fig. 4 



PWTC PWPD PWTD 




,v i . i . i . i . i 

-1,0 -0,8 -0,6 -0,4 -0,2 0,0 

H 





Cluster size 



Cluster size 



